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ABSTRACT 


This study deals with the utilization of 
piezo-electric actuators in controlling the structural 
vibrations of flexible beams. 

A Modified Independent Modal Space Control (MIMSC) 
method is devised to enable the selection of the optimal 
location, control gains and excitation voltage of the 
piezo— electric actuators in a way that would minimize the 
amplitudes of vibrations of beams to which these actuators 
are bonded, as well as the input control energy necessary to 
suppress these vibrations. 

The developed method accounts for the effects that 
the piezo-electric actuators have on changing the elastic and 
inertial properties of the flexible beams. 

Numerical examples are presented to illustrate the 
application of the developed MIMSC method in minimizing the 
structural vibrations of beams of different materials when 
subjected to different loading and end conditions using 
ceramic or polymeric piezo-electric actuators. 

The obtained results emphasize the importance of the 
devised method in designing more realistic active control 
systems for flexible beams, in particular, and large flexible 
structures in general. 



NOMENCLATURE 


b 



D 



E 


i 


f 




F 

r 


8 1 


,2 




K . 


1 


K 


width of beam and piezo-actuator, ■ 

partitioned matrices of the modal shape matrix, i«C or R 
electric charge constant of piezo-actuator, m/v 
distance to neutral axis of composite beam measured 
from its lower edge, m 

Young's modulus of elasticity of piezo-actuator and 

2 

beam respectively, N/m 

Young's modulus of elasticity of element i of the beam, 
N/m 2 

vector of modal forces and moments, N or Nm 

vector of external forces and moments, N or Nm 

modal forces acting on the controlled modes, N or Nm 

physical forces applied to control C modes, N or Nm 

modal force controlling the i^ mode, N or Nm 

vector of external forces and moments acting on element 

i , N or Nm. 

modal forces acting on the residual modes, N or Nm 
physical forces acting on the residual modes, N or Nm 
modal position and velocity feedback gains 
area moment of inertia of actuator and beam about the 

4 

neutral axis of the composite beam respectively, m 
area moment of inertia of the element i, m^ 
mass moment of inertial of the composite beam at node 
l, kg-m 

stiffness matrix of the element i 

overall stiffness matrix of beam-actuator system 


length of element i, ■ 

■ i mass of the composite beam at node 1 , kg 

M mass matrix of beam-actuator system 

M ei external moment acting on i** 1 node of beam, Nm 

piezo-electric moment generated by piezo-film, Nm 
N number of elements of the beam 

R weighting factor 

t^ 2 thickness of piezo-actuator and beam respectively, m 
n. modal displacement of node i, m 
iij modal velocity of node i, m 
U modal coordinates of the flexible system 

v voltage applied across the piezo-electric film, volts 

external force acting on the i^* 1 element of the beam, N 
y^ the linear translation of node i, m 

w, _ beam and film mass per unit length respectively, kg/m 

Oft 

Greek Letters 


f actor=l or zero if actuator is bonded to element i or 
not respectively 

<5. deflection of node i, m or rad. 
x 

6 deflection vector of all nodes of the beam, m or rad. 

2 

<5 acceleration vector of all nodes of the beam, m/s or 

rad/s^ 

piezo-electric strain in piezo-actuator, m/m 
6 i angular deflection of node i, rad. 

X diagonal matrix of the eigenvalues of the system 

o. 0 bending stresses in piezo-actuator and beam 


2 

respectively N/m 

2 

Oj piezo-electric stress in actuator, N/m 

4> modal shape matrix of the eigenvectors of the flexible 

system 

uk resonant frequency at the i t ^ 1 normal mode, rad./s 

Superscripts 

T designates transpose 

-1 designates inverse 



INTRODUCTION 


Active vibration control systems are becoming essential 
and viable means for minimizing the vibrations of large 
flexible structures which are intended to provide stable 
bases for precision pointing in space. Such control systems 
aim at compensating for the inherently low natural damping 
characteristics of the flexible structures in order to 
suppress their oscillatory behavior over a wide band of 
excitation frequencies without compromising their structural 
integrity. In this way, it would be possible to meet the 
strict operational and functional constraints which are 
currently imposed on these structures. 

Distinct among the presently available active control 
systems are those that rely in their operation on 

piezo-electric actuators. Such systems have proven to be 

experimentally effective in controlling the vibrations of 

* 

simple structural elements such as rectangular beams [1-2] 
and hollow cylindrical masts [3]. The effectiveness of these 
systems is coupled also with the light weight, high force and 
low power consumption capabilities of the piezo-electric 
actuators [4-8], These features rendered this class of 
actuators to be an attractive candidate for controlling 
structural vibrations. 


Numbers between brackets designate references at end of 



The present state-of-the-art of this type of actuators 
has been Halted to the analysis and testing of their 
characteristics [9-11] as influenced by their geometrical or 
operational conditions. A recent attempt has been made by Baz 
[12], to select their optimal geometrical parameters and 
location which are best suited for a particular structure 
subjected to known static loading conditions. The developed 
synthesis procedure has proven to be essential to the 
successful integration of the actuators into the structure in 
order to minimize its static deformation. 

However, no effort has been done to optimize the control 
of the vibrations of a multi-mode flexible system using a 
small number of optimally placed piezo-actuators through the 
development of efficient and realistic control algorithm that 
ensures minimal amplitudes of oscillation and input control 
energy . 

It is, therefore, the purpose of this study to devise 
such an optimal control method that is based on a 
modification of the well known Independent Modal Space 
Control(IMSC) method [13-15]. The devised method accounts for 
the spillover from the controlled modes into the uncontrolled 
modes due to the use of fewer actuators than the modeled 
modes. The developed method incorporates also an optimal 
placement procedure that will enable the selection of the 
optimal location of the piezo-electric actuators in the 
structure in order to guarantee a balance between the 



suppression of the amplitudes of vibration and the input 
control energy. 

Furthermore, the developed procedure considers the 
effect that the piezo-electric actuators have on changing the 
elastic and inertial properties of the structure to which 
they are bonded to. Such changes result in modifying the 
normal modes of the structure one way or another depending on 
the location of the actuator. 

In this study, the emphasis will be placed on 
piezo-electric actuator-beam systems to illustrate the 
devised control strategy. But, without any further 
modifications, the developed method can be readily applied to 
large structures. 



THE PIEZO-ELECTRIC ACTUATOR-BEAM SYSTEM 


A. General Layout 

Figure (1) shows a general layout of a flexible beam A 
whose vibration is to be controlled by piezo-electric 
actuator B. The beam, under consideration, can generally be 
made of several steps which are not necessarily of the same 
thickness or the same material* The interfacial nodes between 
the different steps can be subjected to external forces, 
moments or both. Further, the degrees of freedom of any node 
can be limited to linear translations, angular rotations or 
restrained completely depending on the nature of support at 
the node under consideration. 


In this study, the beam is assumed to have rectangular 
cross section of constant width b. The beam is considered to 
deflect in the transverse direction due to the flextural 
action of the external forces and moments. 


In Figure (1), the piezo-electric actuator B is shown 
bonded to the element i of the flexible beam to form a 
composite beam. When an electric field is applied across the 
film, then it will expand if the field is, for example, along 
the polarization axis of the film and will contract if the 
two were out of phase. The expansion or contraction of the 
film relative to the beam, by virtue of the piezo-electric 
effect, creates longitudinal bending stresses in the 


creates 



Piezo-Electric Actuator 



actuator-beam system 




composite beam which tend to bend the beam in a manner very 
similar to a bimetallic thermostat. 

With proper selection, placement and control of the 
actuator, it would be possible to generate enough 
piezo-electric bending stresses to counter balance the effect 
of the exciting forces and moments acting on the beam in a 
way that minimizes its structural vibrations. 

fi. Finite Element Model Of An Actuator-Beam Element 

Figure (2) shows a schematic drawing of a piezo-film A 
bonded to an element B of the flexible beam. 

If a voltage v is applied across the film, a 
piezo-electric strain c £ is introduced in the film and can be 
computed from : 

e £ - (d/t j )*v (1) 

where d is the electric charge constant of the film, m/v 

tj is the thickness of the piezo-electric actuator, m 

This strain results in a longitudinal stress o^ given by : 

o f « EjCd/t^v (2) 

where Ej is the Young’s modulus of elasticity of the film, 




Figure (2) - Schematic drawing of an actuator bonded 

to a flexible beam element 
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N/m 2 

This, in turn generates a bending aoment around the 
neutral axis of the composite beam, given by : 

" (t l +t 2" D) 

M f - / o f (b*y)dy (3) 

-(t 2 -D) 

where t 2 is the thickness of the beam, m 

b is the width of the beam and the piezo-film, m 

In equation (3), D is the distance of the neutral axis from 

the lower edge of the beam which can be determined by 
considering the force balance in the longitudinal direction X 
of the beam, or : 

Jo 1 dA + J° 2 dA = 0 (4) 

film beam 

or 

-(t 2 -D) D 

Ejb J ydy + E 2 b / ydy * 0 (5) 

-( tj+tj-D) -(h 2 -D) 

• 2 

where E 2 is Young's modulus of elasticity of the beam, N/m 

Equation (5) yields the following expression for D : 


Equation (2), (3) and (6) can be combined to determine the 


bending moment generated by the piezo-film on the 

composite beam as follows : 


d*b*(t 1 +t 2 )*(E 1 E 2 t 2 )*v 

2»(E 1 t 1 +E 2 t 2 ) 


(7) 


For this composite beam, it can be easily shown [16] 
that it has a flextural rigidity (E^I^) given by : 


Vi 


E 1 I 1 + E 2 X 2 


( 8 ) 


where 1^ and I 2 are the area moments of inertia of the film 
and the beam about the neutral axis respectively. 


Let us now assume that the composite beam, shown in 

Figure (2), extends a length 1^ between two nodes (i) and 

(i+1). Further, it is assumed that the external forces and 

V. . as well as the external moments M . and M . . are acting 
l+l ei ei+1 ° 

on the beam at nodes i and i+1 respectively. Then, the 
resulting linear and angular deformations of the beam y i and 
6 . as well as and 6 i+1 at the nodes * and i+1 » 

respectively, can be related to the loads acting on the 
element as follows [17] : 
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Equation (9) can be rewritten as : 


Vi 


( 10 ) 


where is the resultant forces and moments vector acting 

on the beam element i, N 

is the stiffness matrix of the composite beam 
element i, N/m 

6 ^ is the deflection vector of the nodes bounding the 
beam element, m 


Equation (9) constitutes the basic finite element model 
that relates the external loads (V and M^) and piezo-electric 
moments (M^) to the deflections (y and e ) of the element as a 
function of its elastic and inertial parameters. 


The equation can be equally used for any element of the 
beam whether it has a piezo-film bonded to it or not. In the 
latter case, is set to zero and flextural rigidity is 
set equal to that of the flexible beam element under 


consideration 
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The force-displacement characteristics of the individual 
elements of the beam-actuator system, as given for element i 
by equation (9), are combined to determine the overall 
stiffness of the beam system. 

The equilibrium conditions of the overall structure will 
be expressed as : 

External forces and moment _ forces and moments 
acting on the nodes = acting on the elements 

of the overall system at these nodes 

or 

N+l N+l 

F ■ X r i - A Vi * « < n > 

1=1 1=1 

where K is the overall stiff nessmatrix of the system (2n*2n) 

Bathe and Wilson [18], Yang [19] and Fenner [20], for 
example, show how to generate the overall matrix I from the 
stiffness matrices of the individual elements. 

The inertial properties of the composite actuator-beam 
system are determined using the lumped mass method where the 
mass and rotational inertial of each element is distributed 
among the nodes bounding the element [17]. 


Therefore, the diagonal mass matrix M (2n*2n) for the 
actuator-beam system, shown in Figure(l), can be written as : 
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where « K w b +Y l w f )* L j ]/2 

J 1 ■ I<W£>* L 1 3 ]/ 2 

"i ' f( V Y i-iV* L i-i +< VVf >#L i ]/2 
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'b,f 


beam and film mass per unit length 
respectively, kg/m. 

1 if an actuator is bonded to beam element i< 
0 if not. 


The stiffness and mass matrices K and M, defined by 
equations 11 and 12, are used to define the dynamic equations 
of motion of the actuator-beam system. 



The equations of motion of the actuator-beam system can 
be written as follows : 

M2 + - F (13) 

where <$ is the acceleration of the nodal points of the 

2 

composite beam, m/s 

Equation (13) is put in the self-adjoint modal space by 
using the following weighted modal transformation [21] : 

6 « *U (14) 

where U is the modal coordinates of the system 

4> is the weighted modal shape matrix of the 
eigenvectors of the flexible system 

Using such transformation, reduces the coupled equation 
of motion (13) to the following uncoupled form : 

U + XU = f (15) 

where x is a diagonal matrix of the eigenvalues of the system 


f is the modal force matrix given by 



The nodal representation of the flexible bean-actuators 
system la utilised, as a basis for this study, as it nakes 
the design and the selection of the gains of the active 
control systems rather sinple as compared to the coupled 
system representation. This is due to the fact that the 
design of modal controllers using, for example, the 
Independent Modal Space Control (IMSC) method [13-15] has 
been proven, in numerous studies such as [15], to be 
computationally more efficient than the conventional Riccati 
or Pole Assignment methods. 


CONTROL STRATEGY OF ACTUATOR-BEAM SYSTEM 

The control strategy used in controlling the vibration 
of the flexible beam system is based on a modified version 
(MIMSC) of the Independent Modal Space Control (IMSC) method. 
Details of the method are presented in Reference [22]. In 
essence, the MIMSC utilies small number of optimally placed 
actuators to control a large system that has large number of 
degrees of freedom. The method accounts also for the effect 
of interaction between the controlled and uncontrolled modes. 

Considering the modal representation of the elastic 


system, one could rewrite the modal forces f as : 
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W 
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(17) 


where ^ are the modal forces on the controlled and 

residual modes respectively • 

£ are the physical forces on the controlled and 
residual modes. 

^(lj) is the modal shape at mode i and location 1^. 


The above equation can be rewritten as : 


B cc b cr1 | f c 


b rc b rrJ L f r 


(18) 


If only C modes are controlled with equal number of control 
forces F^ f then F^eO and equation (18) reduces to : 


f s B F 
E C D CC C 


( 19 ) 


P c R p 

r R D RC'c 


In the IMSC method, it is assumed that the control 
forces F c will not contribute to the excitation of the 
residual higher order nodes* Accordingly* it was assuaed that 

This of course can only be true if the nunber of 
controlled nodes is very large compared to the number of 
residual modes or when the residual modes are at much higher 
frequency band than the controlled modes. If these two 
conditions are not satisfied* then there will be considerable 
interaction between the controlled and residual modes. 

The MIMSC method considers such interaction by 
calculating the optimal modal control forces [f^,] using the 
IMSC close form solution of the Riccati Equation such that 
the control force f^ of the i t * 1 mode, as given by [14], is : 

f i * + 8 2 “i^ /R (21) 

where R is a factor that weighs the importance of 

minimizing the vibration with respect to the 
control forces. 

'* ) i is the resonant frequency at the i tfl normal 

mode . 

u.^, are the modal displacement and velocity 

respectively. 

gj, g 2 are the modal position and velocity feedback 
gains given by [22] as : 
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(23) 


Accordingly, the displacement Uj and Telocity a^ at the 
i** 1 aode can be feedback and used along with equations (21) t 
(22) and (23) to determine the modal control force f^. 


Once these forces are calculated, equation (19) is 
solved to give the physically applied control forces Fq as 
follows : 



(24) 


Then equation (20) is used to calculate the modal forces 
fg that would excite the residual modes and which are 
generated by the spillover from the controlled modes. 
Definitely these fg are not equal to zero as originally 
assumed in the IMSC method. 


Equations (17) and (15) can then be integrated with 
respect to the time to determine the modal displacements (u^ 
and velocities (u^) which can, in turn, be used again to 
compute the modal forces f and so on. 

From the modal displacements and velocities, the 
physical state ( 4) of the flexible system can be determined 
from equation (14). A relationship can therefore be 
established between the physical state 6 of the system and 
the physical control forces F^ applied to it. 


It is very important to point ont here that the 
magnitude of the nodal forces fg depends primarily on the 
magnitude “^*s of the controlled modes as veil as the modal 
state variables u and u. On the other hand, the magnitude of 
the actual physical control forces Fg depends mainly, for a 
given controlled mode, on the point of application of these 
forces as defined by the matrix Bgg~*. Therefore, minimizing 
fg does not necessarily means that Fg vill be minimum lnspite 
of the fact that it is represented as a linear combination of 
fg. This is simply because the coefficients of the linear 
combination, which are elements of the Bgg“* matrix, depend 
on the placement strategy of the control forces Fg. One 
could still find an optimally placed set of physical control 
forces Fg such that the physical displacements and control 
forces would assume minimum value. 

This optimum placement of the physical control forces is 
an important feature of the MIHSC method and vill be 
demonstrated to be essential part of the design of the active 
control system. 

One should stress here also that if all the modes are 
controlled and there is no residual modes then the conditions 
for minimizing fg will make Fg minimum as well. But, in real 
large structures this vill be unlikely to happen as the 
number of controlled modes is much smaller than the number of 
modeled modes. Therefore, it is essential to augment the IMSC 
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■ethod with an optimal placeaent algorithm to guarantee 
efficient design of the control system. 

The optimum placement of the actuators is implemented 
through the use of the Univariate Search method [23] which 
varies the location of one actuator at a time in order to : 

Minimize / (& 2 + RF C 2 ) dt (25) 

In other words, the optimal placement algorithm 
minimizes the weighted sum of the amplitudes of vibration and 
the generated control forces. The weighting factor R is 
selected by the designer to emphasize the importance of 
damping out the vibration over the expended control energy 
(when R<<1) or vice versa when R>>1 . Equal importance of the 
two parameters is achieved with R=1 . 

The MIMSC method incorporates also an extremely 
important feature which is based on the "TIME SHARING" of a 
small number of actuators in the modal space to control large 
number of modes. Details of such important feature are 
presented in References [29-30], In few words, the Time 
Sharing Control Strategy generates modal control forces to 
control, at the first time interval, the first through the 
C*"* 1 modes using C actuators then, at the second time 
interval, it provides signals to control the second through 
the (C+l) t ^ 1 modes followed by commands to control the third 
through the (C+2) t * 1 modes and so on until all the modeled 


■odes are controlled in this sequential fashion. Once all the 
modeled modes have received their share from the control 
action the cycle is repeated again to effectively damp out 
all the modes of vibration with few number of actuators. This 
strategy has been shown to result in efficient control of the 
vibration of large structures with relatively small number of 
actuators when other existing methods fail to do so [22]. 

Figure (3) outlines a flowchart of the MIMSC method 
indicating the main steps of optimal placement and time 
sharing of the actuators as well as the consideration of the 
spillover between the controlled and residual modes. 


Application Of MIMSC Method 
To Piezo Actuators-Beam System 

The MIMSC method is utilized to design active vibration 
controllers for flexible beams of particular configuration 
when subjected to specific external loading and end 
conditions. The controllers rely in their operation on one or 
more piezo-electric actuators which are optimally selected, 
placed and controlled without violating their polarization or 
strength constraints. Two types of piezo-electric actuators 
are considered in this study namely : ceramics and polymeric 
actuators whose physical properties are given in Table (1). 




Figure (-3) - Flov chert of the MXMSC Computational Algorithm 



Table (1) - Properties of piezo-electric actuators 


Actuator Type 

Ceramlc(Ref .8&10) 

Polymeric (Ref .24) 

Material 

PZT 

G1278 

KYNAR (PVF2) 

Charge coefficient d 

123 

270 

23 

10 * 2 (m/v) 




Young's modulus 
(G N/m 2 ) 

139 

60 

2 

9 

Max. voltage v 
(Mv/m) 

1 

0.7 

30 

Max. tensile 
strength 

45 

47.5 

33 - 55 

(M N/m 2 ) 




Density 

(kg/m 3 ) 

7500 

7400 

1780 


NUMERICAL EXAMPLES 

The application of the MIMSC aethod is illustrated by 
considering a flexible straight beam, shown in Figure (4), 
which is 0.0125m wide, 0.0021m thick and 0.15m long. The beam 
is modeled by three-element model for the sake of 
demonstrating the MIMSC method. 

The effect of varying beam loading and end conditions as 
well as beam and actuator materials on the optimal 
configuration of the system is considered in great details. 
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A. Vibration of Cantilever Beam 

The flexible bean, under consideration, is restrained at 
node 1 to form a cantilever beam whose free end, node 4, is 
subjected to either transverse impulsive loading or 
sinusoidal excitations. The considered impulsive load is 
selected to have 8 magnitude of 0.1N and duration of 1.0ms 
whereas the sinusoidal excitations are assumed to have the 
same magnitude of 0.1N and a frequency ranging between 15 to 
1000 Hz. 

The effect of the two types of loading on the time 
and/or the frequency responses of the piezo electric 
controlled beam is determined using the MIMSC method and 
compared to the corresponding characteristics of the 
uncontrolled beam. Such comparison are used to determine the 
effectiveness of piezo-electric actuators in damping out the 
vibrations of the undamped flexible beam. 

(I). Response to The Impulsive Loading 

Figures (5-a) and (5-b) show the time histories of the 
amplitudes of transverse vibrations of the uncontrolled beam 
when it is made of steel and aluminum respectively. It can be 
seen that the beams will continue to vibrate in a limit cycle 
fashion even after termination of the input impulse. This is 
due to the fact that the beams are assumed to have no 
inherent structural damping. 




The figures indicate also that the alunlnun bean will 
experience higher amplitudes of vibrations than the steel 
beam . 


The amplitudes of vibrations displayed in Figures (5-a) 
and (5-b) are used as a datum for measuring the effect of the 
thickness, material, location and number of piezo-electric 
actuators on controlling the vibration of the cantilever beam 
as will be described in that follows. 

(a). Effect of Actuator Thickness 


Figures (6-a), ( 6— b ) and (6-c) illustrate the time 
histories of the amplitudes of transverse vibrations of the 
steel cantilever beam when controlled by the MIMSC method 
using one piezo-electric actuator that has a thickness of 
0.000525m, 0.00105m and 0.0021m respectively. In 
dimensionless form, the actuators are selected to have 
thicknesses equal to 1/4, 1/2 and 1 that of the beam 
thickness respectively. All the considered actuators are made 
of PZT4 ceramic are placed between nodes 1 and 2, i.e. bonded 
to the element near the fixed end of the beam. 

Figures (7-a), ( 7— b ) and (7-c) display the corresponding 
control voltages generated, according to the MIMSC method, to 
power the considered four actuators respectively. 







Figures (6-a) through (6-c) indicate that using one 


actuator which is capable only of producing bending Moments, 
it was possible to control all the three linear translations 
and the three angular rotations of the cantilever beam. This 
feature is an important feature of the MIMSC method where one 
actuator is shared to control more than one mode of vibration 
of the flexible structure. 

Detailed analysis of the figures is summarized in Table 
(2) to Indicate the effect of the thinkness of the 

piezo-electric actuator on the displacement index U^, control 
force index U and control energy index U v . These Indices are 
defined by the following expressions: 

t=t # N , 

U. = 2 2 6 . * (26) 

t-0 i=l 

t=t* N 0 

U = 2 I F/**t (27) 

t«0 i*l 

t=t # N 

Up =2 2 |«i*Fj At (28) 

t=0 i=l 

where At is time increment for integrating the equations of 
motion taken as : 


At « 1/10*2*3. 1416/wmax (29) 

to maximum natural mode of vibration of the elastic 

max 


t maximum time of integration, sec 
N number of degrees of freedom of the structure 

Table (2) - Effect of actuator thickness on the displacement, 
control force and control energy indices for a 
steel cantilever beam controlled with one PZT 
actuator at element 1 


Thickness (mm) 


0.525 

1.050 

2.100 

4.200 

8.400 

Displacement Index, 

xlO 9 

2.264 

1.750 

1.354 

1.195 

0.975 

Control force Index, 

xlO 5 

4.740 

5.299 

6.003 

10.337 

3616.300 

Control Energy Index 

, xlO 6 

0.354 

0.343 

0.350 

0.307 

28.690 

Table (2) 

indicates 

clearly 

that 

incr 

easing 


thickness of the actuator results in decreasing the amplitude 
of vibrations of the beam considerably as measured by the 
displacement index U^. Such a reduction is attributed 
primarily to the stiffness effect that result from bonding 
thicker actuators to the original flexible beam. However, 
using such thicker actuators would require large control 
forces to damp out the vibrations of the beam-actuator 
system. This is reflected directly into high control voltages 
as can be seen clearly from Figures (7-a) through (7-c). 

Therefore, we have a push and a pull situation where 
thicker actuators are preferred to minimize the amplitudes of 
oscillation and thinner actuators are desirable if one is to 
limit the magnitude of the control forces and control 
voltages. Accordingly, a balance can be achieved between 


minimizing the displacement index without having to produce 
excessively large control forces. This balance can be 
attained, if we consider the control energy index Ug which is 
a measure of the energy spent to control the vibration of the 
beam. From Table (2), it can be seen that there is an optimum 
actuator thickness, t=4.2mm, at which the control energy 
index assumes a minimum value. 

One should also point out here that using the 
piezoactuators results in damping out completely the 
vibrations of the beam when compared to the case of the 
uncontrolled beam as can be seen from a comparison between 
Figures (6-a), (6-b), (6-c) and (5-a) respectively. 


It is 

important 

to 

note here 

that bonding 

the 

different 

thickness 

actuators 

to 

the beam 

is found 

to 

result in 

changing 

the natural 

modes of 

vibration 

of 

the beam 


significantly. For example, with an actuator thickness equals 
to 1/4 the beam thickness the first normal mode occurs at 
67Hz but when the thickness is increased to the beam 
thickness, the first mode decreased to 55.3Hz. Therefore, It 
is essential to consider the effect of introducing the 
actuators on changing the elastic and inertial properties of 
the flexible system before designing the active control 
system. 

Also, it is essential to note that the design of the 
active modal controller is carried out with a weighting 


factor R«50 in order to maintain the control forces low 
enough to avoid depolarization but high enough to control 
effectively the beam system. Detailed account for the effect 
of the weighting factor R on the dynamics of the beam system 
will be given later on. 

(b) Effect of Number of Actuators 

The effect of the actuator thickness should also be 
considered in view of the number of actuators to be used in 
controlling the vibration of the steel beam. For example, 
considering the case of two actuators to be used to control 
the six degrees of freedom of the beam. The considered two 
actuators are made of PZT of the same thickness and placed at 
elements 1 and 3. The effect of varying the thickness of the 
two actuators from 1/4 to 4 times the beam thickness on the 
vibration of the beam is shown in Figures (8-a) through (8-c) 
respectively. The corresponding control voltages are 
displayed in Figures (9-a) through (9-c). 

Figures (8-a), (8-b) and (8-c) indicate that increasing 
the thickness of the actuators produced favorable damping 
characteristics as it results in reducing the maximum 
amplitude of transverse vibration of the beam. But this is on 
the expense of a significant increase in the required control 
voltages. For example, increasing the actuator thickness from 
1/4 to 4 times the beam thickness reduces the maximum 
amplitude of vibration from 2.05E-5m to 1.06E-5m whereas it 
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increases the maximum control voltage from 5.5v to 720v. 
Therefore, the benefits gained bj the reduction in the 
amplitude of vibration are by far less than the drawbacks 
resulting from excessive control voltages. A proper balance 
between these two parameters is essential for an optimally 
controlled system. Such a balance can be attained by 

considering the control energy index listed in Table (3). 

Table (3) - Effect of actuator thickness on the displacement, 
control force and control energy indices for a 
steel cantilever beam controlled with two PZT 
actuators at elements 1 and 3 


Thickness (mm) 

0.525 

2.100 

4.200 

8.400 

§ 

Displacement Index, xlO 

0.628 

0.498 

0.412 

0.326 

0 

Control force Index, xlO 

12.140 

10.000 

33.250 

5346.800 

Control Energy Index, xlO^ 

0.675 

0.515 

1.409 

70.900 


Table(3) suggests that using two actuators each of a 
thickness equal to that of the original beam results in 
minimal control energy expenditure. 

Furthermore, comparing the results listed in 
Tables(2)and(3) one could see that it is more beneficial, 
froc the displacement and control energy indices, to use the 
two optimal actuators to control the beam vibration rather 
than using just one optimal actuator. This is inspite of the 
fact the total weight of the two optimal actuators is equal 
to the optimal single actuator. 


Therefore, for a given weight of the structure and the 
piezoactuators, it is advantageous to distribute the 
piezo-material over the structural element to minimize its 
vibrations . 

(c) Effect of Actuator Location 
(i) Single Actuator 

Figures (10-a), (10-b) and (10-c) show the time 
history of the amplitude of transverse vibration of the steel 
beam when a PZT piezo-electric actuator, of a thickness twice 
the beam thickness, is placed at element 1, 2 or 3 

respectively. The corresponding control voltages are shown in 
Figures (11-a), (11-b) and (11-c) for the considered three 

placement strategies respectively. 

The associated displacement, control force and control 

energy indices are listed in Table (4). 

Table (4) - Effect of actuator location on the displacement, 
control force and control energy indices for a 
steel cantilever beam controlled with one PZT 
actuator of thickness=4 . 2mm 


Actuator Location 

Element 1 

Element 2 

Element 3 

0 

Displacement Index, xlO 

1.195 

1.516 

1.326 

0 

Control force Index, xlO 

10.337 

356.75 

2658.800 

Control Energy Index, xlO^ 

0.307 

11.07 

61.800 
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It is i therefore, obvious that when a single actuator is 


to be used to control the vibrations of the cantilever beam, 
than this actuator should be bonded to element 1 near the 
fixed end of the beam. Such a placement strategy is optimal, 
from all points of views, as it minimizes the displacement, 
control force and control energy indices all at the same 
time . 


(ii) Two Actuators 


When two PZT actuators are used to control the steel 
cantilever beam, the displacement, control force and control 
energy indices are as listed in Table (5). 


Table (5) - Effect of actuator location on the displacement, 
control force and control energy indices for a 
steel cantilever beam controlled with two PZT 
actuators of thickness=4 . 2mm 


Actuator Location 

Element 
1 and 2 

Element 
1 and 3 

Element 
2 and 3 

5 

Displacement Index, xlO 

0.342 

0.412 

0.393 

9 

Control force Index, xlO 

27.170 

33.250 

678.360 

Control Energy Index, xlO^ 

2.100 

1.691 

31.170 


The obtained results indicated that it is preferable to 
place the actuators at element 1 and 3 in order to obtain the 
minimum displacement, control force and control energy 



(d) Effect of Weighting Factor R 


(1) Single Actuator 

Figures (12-a), ( 1 2— b ) and (12-c) show the time history 
of the amplitude of transverse vibration of the steel 
cantilever beam, controlled by one PZT actuator placed at 
element 1, for values of the weighting factor R of 1, 20 and 
35 respectively. The associated time histories of the control 
voltages are shown in Figures (13-a), ( 1 3— b ) and (13-c) 
respectively. A summary of the displacement, control force 
and control energy indices is given is Table 6 for values of 
the weighting factor between 1 and 200. 


Table (6) - Effect of the weighting factor R on the 

displacement, control force and control energy 
indices for a steel cantilever beam controlled 
with one PZT actuator of thickness»4 ,2mm placed 
at element 1 


Weighting Factor R 

1 

20 

35 

50 

100 

200 

G 

Displacement Index, xlO 

0.204 

0.778 

1.017 

1.195 

158.26 

199.20 

_ q 

Control force Index, xlO 

2450 

25.200 

14.100 

10.300 

6.18 

4.027 

Control Energy Index, xlO^ 

9.017 

0.497 

0.386 

0.308 

0.227 

0.173 


The presented figures and the obtained results indicate 
that emphasizing the importance of the control force over the 
amplitude of vibration by increasing the value of the 
weighting factor R results in lowering the magnitude of the 
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control force and slowing accordingly the process of damping 
out the structural vibrations. For example, placing equal 
weigh on the displacement and control force results in 
maximum amplitude of displacement of 1.26E-5m, as seen from 
Figure (12-a), but this would require a maximum control 
voltage of 236v. Increasing the weighting factor R to 35 
results in a slight increase in the maximum amplitude of 
vibration to 1.47E-5m but this is accompanied with a 
considerable drop in the maximum control voltage to 11. 3v as 
can be seen from Figures (12-c) and (13-c) respectively. 
According, increasing R, up to 50, results in significant 
reduction in the control voltage without producing serious 
degradation in the damping out of the beam vibration. Beyond 
R®50, it is found that vibration damping process is becoming 
excessively slow to the extent that overshadows the savings 
in the control force or the control energy indices. 

(ii) Two Actuators 

When two PZT piezo-electric actuators are placed at 
elements 1 and 3 of the steel beam, the resulting time 
history of the amplitudes of transverse vibrations are shown 
in Figures (14-a) and (14-b) for values of the weighting 
factor R of 1 and 20 respectively. Figures (15-a) and (15-b) 
show the corresponding control voltages. 

A summary of the effect of R on the displacement, 
control force and control energy indices are given in Table 
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Table (7)- Effect of the weighting factor R on the 

displacement, control force and control energy 
Indices for a steel cantilever beam controlled by 
two PZT actuators of thickness=4.2mm placed at 
elements 1 and 3 


Weighting Factor R 

1 

20 

50 

0 

Displacement Index, xlO 

0.0504 

0.259 

0.4125 

0 

Control force Index, xlO 

23948 

119.450 

33.2500 

Control Energy Index, xlO^ 

29.520 

2.253 

1.4090 


Table (7) indicates that increasing R results in 
reducing the control force and the control energy but at the 
same time produces higher values of the displacement index. 
Again, as mentioned before, using two actuators is seen to be 
more effective in damping out the beam vibrations. 

(e) Effect of Beam Material 

The effect of changing the material of the cantilever 
beam, under consideration, from steel to aluminum on the 
resulting amplitudes of transverse vibrations and the voltage 
required to control it with one PZT actuator is shown in 
Figures (16-a) and (16-b) respectively. The actuator has a 
thickness of 4.2mm and is placed at element 1. The controller 
is designed with a weighting factor R=50. For this case, the 
displacement, control force and control energy indices are 
computed to be 10.4E-9, 19.85E-9 and 0.928E-6 respectively as 


CSIBGI^AL PAGE 13 

Qf POOR QUALITY 


4.XMS 




in 


I.K«« I.W-K 1.4KB l.lil-K l.«Ml 

HU - src. 


vib ratio© of as ©antllcrer hn tOBtoUtd tr otm pb 

ae tutor plaoU at aluut 1. V — WB 


4.31f‘£l 


(‘i'JUMI li.“ ’( ;M> -• 

rc tcicMi'Fs: n:ns 1 

f(i):8.1H WCT.rfc. 

toe 


2.47M1 


-l.nuti 



ImLLUAiu^ 


TIKE - sec. 


ruur* (16-b) - rime history or the ©octroi aclt ft raouirad 
to control an aliminus eaxstilarer b*u by ooe FIT actuator pUcwd 
at elaatfit 1. 


compared to 1.19E-9, 10.34E-9 and 0.307E-6 for the case of 
the steel beam. Therefore, it is evident that the selected 
PZT actuator is more effective in damping out the vibration 
of the steel beam than the aluminum beam. 

(f) Effect of Actuator Material 


The effect of using a Kynar polymeric actuator, instead 
of PZT actuator, on the amplitudes of vibration and the 
voltage required to control the vibration of the aluminum 
cantilever beam is show in Figures (17-a) and (17-b) 
respectively. The Kynar actuator is assumed to have a 
thickness of 4.2mm and be placed at element 1. 

A comparison between Figures (16-a) and (17-b) indicates 
that the maximum amplitude of vibration has increased from 
4.36E-5m to 4.88E-5m when the material of the piezo-actuator 
is changed from PZT to Kynar. This amplitude increase is 
about 11.9%. But, the required maximum control voltage is 
found to increase from 63.1 to 141 volts respectively. This 
amounts to an increase of 123.4%. Therefore, from both the 
displacement and control voltage, it is preferable for the 
arrangement under consideration to use PZT actuators. 

A detailed analysis of the obtained results indicate 
that the displacement, control force and control energy 
indices are 15.28E-9, 5.86E-9 and 0.964E-6 respectively. 



II. Response to Sinusoidal Loading 


The steel cantilever bean, which is controlled bj one 
PZT actuator of thickness 4.2mm placed at element 1, is 

subjected to external sinusoidal transverse force acting at 
node 4. The applied force has a magnitude of 0.1N and 
frequency ranging from 15Hz to 1000Hz. 

Figures (18-a), (18-b) and (18-c) show a comparison 
between the time histories of the amplitudes of vibration of 
the controlled and uncontrolled beam at frequencies of 100, 
900 and 3100 rad/sec respectively. 

From these figures, it can be seen that the 
piezoactuator is capable of attenuating the response of the 
beam to the sinusoidal excitation throughout the considered 
range of frequency. The attenuation is however more 
significant at high frequency than at lower frequencies. 

A summary of the effect of the excitation frequency on 
the ratio of the amplitudes of vibration of the controlled 
beam to those of the uncontrolled beam, in decibels, is shown 
in Figure (19). Figure (19-a) displays the linear amplitude 
ratios and Figure (19-b) shows the angular amplitude ratios 
for the different nodes. 


B. Vibration of Overhung Beam 
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The three element steel beam, under consideration, is 
fixed at node 1 and is restrained from linear transverse 
motion at node 3. Such configuration is considered as an 
example to illustrate the effect that the end conditions of a 
flexible structural element have on its vibration control 
character 1st ices. 

The effectiveness of the control strategy will be 
compared to the time history of the amplitudes of vibration 
of the uncontrolled beam which are shown in Figure (20). 

The considered overhung beam is assumed to be controlled 
by two PZT actuators each has a thickness of 4.2mm. The 
effect of placing the two actuators at elements 1 and 2, 1 
and 3 or 2 and 3 on the time history of the amplitudes of 
transverse vibration is shown in Figures (21-a), (21-b) and 
(21-c) respectively when the beam is subjected to an impulse 
of 0.1N for a duration of 1ms at node 4. The corresponding 
control voltages are shown in Figures (22-a), (22-b) and 
(22-c) respectively. 

The associated displacement, control force and control 
energy indices are listed in Table (8). 




im-tt 


nmt* r»* tmn' 

gonus ■ 


I.M-K 


-J.W-B 


l.tfttt M«-K t. 32 E -*2 I.Cl-C 
IlhC • !K. 

ri*nra (?i-a) - rw tutor? or to* aaputuu* «r tiumm 

Titration of as orarbunc atari kw abon oontrollnd 07 too m 
actuator* of thloknaana0.0M2 a plaoad at alaaasta 1 aad X, 


1.1IM6 


un-a 


l.fct«M :.W-a 5.27K2 7.JU-K USE-6! 

UK - Sat. 

flfura (21-0) - Tlaa hlatorj of tka aaplltada* of tranararaa 

alt rat lot of an ovtrbuae atoel boaa abac oootrollae toj tao f2T 
aetuatora of tniokna*asG.00t2 a plaoad at aUBanta 1 aad J. 




MUSIC H0> ttlEEL' 
fa idlaKi (I Effitt 02 
WKls.t. 

IriS 


hU .* , 


OF PO <M QUAUTT 


8.UM8 2.7IK2 5.tfH2 I.Jir-K i.UF-Ci 

l!K • stc. 


Tigorm (21«<t) - tim Idatory of tki m^UUAm af t nairm 
ttbratloa of an ovartaon | at Ml toaa vfcta aoctroUad toy tao fZT 
aetuatora of tbloknaa»«0.0042 a plaoad at alaaafits 2 aad 3 . 



ns - sk. 

Fi*ur» (22-t) - Ti»* Mitory of tte ooetrol volU|M required 
to ©octroi u ororbunt otosl Mm by too «T cetoetors of 
tbick&s*»*0.00*2 m ploood at oloooou Icdt. 


5.1CMH - 

nnc fltDKBNC EM (SHE) 

ra tcmm it EUKn is 
HM.W M.Wjk. 

1.MM1 " 

Ma .... 

r 

•MS* - 

r 

n 

C i 
**» • 

m 

• 

1 

f.CK 

■« i.K-n 5.nr-rr mu-k i.rmi 


nr *. i»t. 

Flfurt (2?-b) * Tie* bistory of tbc ooetrol Toltafas required 

to cootrol an otcrtouaf otosl Mss by too PZT oc too tors of 

tM nir nos ■ *0.0042 ■ Heel at oloeseU 1 «ed 3. 

I.2SI.K - 

onec E“ (JTHD 

F37 iCTItHCf :T DiMNT 
FW4.1H W.ttlst:. 

! t54 

1.UI4M - 




M' 

-3.RDN - 

-7.71I+K j 


2.7II-K 5.55 £-62 I.S3E-E l.UI-fi 


IW - MC. 


ri*uro (22-o) - TIm tlotorj of tta ooetrol velt^os rognlrod 
to ooetrol on ortrteoc etool Mm by too PZT actuators of 

t±IckOM.«0.00*2 > pioood at olsents 2 mA 3. 


Table (8) - Effect of placement strategy of two PZT actuators 
of thickness>4 .2am on the displacement , control 
force and control energy indices of an overhung 
ateel beaa 


Actuator Location 

Element 
1 and 2 

Element 
1 and 3 

Element 
2 and 3 

6 

Displacement Index, xlO 

0.0219 

0.0226 

0.0221 

b 

Control force Index, xlO 

7.2520 

0.6190 

0.1664 

Control Energy Index, xlO^ 

9.8900 

1.9870 

0.9140 

The results listed in 

Table 8 

suggest that 

placing the 


actuators at elements 2 and 3 results in minimizing the 
control force and control energy indices simultaneously. 


Conclusions 


This study has presented a comprehensive analysis of the 
control of vibration of simple structural elements using 

piezo-electric actuators. The developed vibration control 
algorithm is based on a Modified Independent Modal Space 
Control (MIMSC) method which accounts for the interaction 
between the controlled and uncontrolled modes, optimal 
placement of the actuators as veil as time sharing of few 
actuators between the modes of the structure. The method is 
found to be extremely effective in controlling all the modes 
of vibration of simple structural elements with very small 
numbers of actuators. 

The MIMSC is developed to account also for the effects 
that the piezo-electric-based active control system it is 
essential to determine the optimal actuator geometry, control 
voltage and location for a particular loading and end 
conditions. It was concluded that there is an optimal 
actuator thickness of which the displacement, control force 


and control 

energy can 

be 

minimized. Thin 

actuators are 

ineffective 

in damping 

out 

the vibration 

whereas thick 


actuators would require much higher control voltages. 

It is observed also that spreading the piezo-electric 
material over the structural element results in lowering the 
amplitude of vibration considerably. Furthermore, the study 
indicated also the importance of selecting the proper 



weighting factor between the damping of vibration and the 
expenditure of control energy. 

Using ceramic PZT actuator is seen to be more effective 
than Kynar actuators in controlling the vibrations of the 
considered simple structural elements without requiring 
excessively high control voltages. 

The study has demonstrated also the effectiveness of the 
developed MINSC method in controlling the vibration of 
flexible beams made of different materials and subjected to 
different types of end constraints. 

The presented results emphasize the importance of the 
MIMSC method in designing more realistic active control 
systems for flexible beam in particular and large flexible 
structures in general. 

The extension of the developed algorithm to large 
structures can be readily done without any modifications. 
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